Theory of Electronic Raman Scattering in Nearly Antiferromagnetic Fermi Liquids 
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A theory of electronic Raman scattering in nearly antiferromagnetic Fermi liquids is constructed 
using the phenomenological electron-electron interaction introduced by Millis, Monien, and Pines. 
The role of "hot spots" and their resulting signatures in the channel dependent Raman spectra 
is highlighted, and different scaling regimes are addressed. The effects of self consistency and 
vertex corrections are compared to a perturbative treatment of the interaction. The theory is then 
compared to Raman spectra taken in the normal state of cuprate superconductors, and it is shown 
that many features of the symmetry dependent spectra can be explained by the theory. Remaining 
questions are addressed. 

PACS numbers: 74.25. Jb, 71.27.+a, 78.30-j 



I 

o 
o 



> 
as 
cn 
o 



o 
o 



X 



I. INTRODUCTION. 

Recent results from angle resolved photo-emission 
(ARPES) studies have shown that the planar quasipar- 
ticlc (qp) dynamics in the normal state of high temper- 
ature superconductors are extremely anisotropic [0 . For 
optimally doped samples with the highest Tc's, the qp 
spectral function observed for Fermi surface crossings 
along the (tt, 0) — (tt, tt) direction in the Brillouin zone 
(BZ) is smeared by strong scattering compared to the ob- 
served spectrum along the (0, 0) — (tt, tt) direction. This 
anisotropy becomes more pronounced as the materials 
are further underdoped with concomitantly lower tran- 
sition temperatures into the superconducting state 
These effects have been ascribed to a strongly momen- 
tum dependent qp self energy E(k, w) whose imaginary 
part is largest for k near the BZ axes and yields a broadly 
incoherent spectral function. Understanding the nature 
and origin of this anisotropy could lead to a better un- 
derstanding of the complex dynamics which cause the 
normal state of the cuprates to be very far away from 
being "normal". 

Many theoretical efforts have been conducted to un- 
derstand this planar anisotropy by calculating S(k,w) 
for various types of interactions and via several meth- 
ods. These efforts usually have focused on understand- 
ing the results from photo-emission. Since ARPES is a 
single particle probe, the spectral function can be com- 
pared directly with experiment and provides direct con- 
tact to the self energy around the BZ. However, most 
other probes (e.g. penetration depth, conductivity, tun- 
neling density of states, etc.) involve averages around 
the BZ and thus are not so sensitive to the effects of 
anisotropic self energies. Electronic Raman scattering is 
an exception. As demonstrated in the superconducting 
state 1^, electronic Raman scattering is a two particle 
probe which has proven to be a useful tool to understand 
qp dynamics on selected regions of the BZ. The matrix 
elements for Raman scattering can be tuned by orient- 
ing incoming and outgoing photon polarizations, thereby 



projecting out different regions of the Fermi surface, fn 
Ref. it was shown how this tuning could be used to 
provide a detailed probe of the momentum dependence 
of the superconducting energy gap. 

In this regard, Raman scattering can be considered 
a complementary probe to ARPES. Raman measure- 
ments on several systems have shown a remarkable chan- 
nel dependence in the normal state for the optimally- 
and under-doped cuprates in a variety of systems (Q-HI- 
For scattering geometries of Big symmetry, an almost |§] 
temperature independent Raman spectrum has been ob- 
served up to frequencies of a few J, while for i?2g ori- 
entations the spectra show a temperature dependence 
that tracks that of the DC resistivity Moreover, 
the overall magnitude of the electronic continuum is ob- 
served to be relatively independent of doping for the i?2g 
channel, but decreases dramatically with undcrdoping for 
the Big channel, with spectral weight transferring out to 
two-magnon energies [0,^. One is to conclude that the 
qps near (tt, 0) which are probed by the Big channel are 
particularly incoherent and become more incoherent with 
underdoping, in agreement with the results from ARPES. 
Coupled to the flatness of the bands there, these qps are 
likely not to participate in transport |^ . Currently there 
is no theoretical understanding of Raman scattering in 
the normal state of the cuprates, as all previous theo- 
ries failed to predict large differences between these two 
channels at low frequencies. This is due to the neglect of 
strongly anisotropic electron interactions. It is the pur- 
pose of this paper to address these issues. 

It is widely believed that strong antiferromagnetic 
(AF) correlations are an important ingredient needed 
to describe the unusual properties found in the normal 
state of the cuprate superconductors and provides the ori- 
gin of the strongly anisotropic dynamics in the cuprates. 
Large qp scattering for momentum transfers near the AF 
wave-vector Q = 7r/a,7r/a have been invoked to explain 
the photo-emission results [^,0. Therefore we have 
studied the phenomenological Nearly Antiferromagnetic 
Fermi Liquid (NAFL) model |l2[ as an example of highly 
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anisotropic scattering to investigate the two-particle Ra- 
man response. The NAFL model has at its center a 
strong magnetic interaction between planar qps which 
is peaked at or near Q and yields anisotropic self ener- 
gies. This model has been wide ly app lied to investigate 
both spin and charge properties . Our results show 

that the Raman response is an ideal probe to clarify the 
nature of qp dynamics on regions of the Fermi surface. 

The outline of the paper is as follows: Section II re- 
views and outlines a generalized approach to calculate 
non-resonant Raman scattering in correlated electron 
systems and introduces the model interaction. Section 
III investigates the NAFL model using simple pertur- 
bation theory. The scattering rate for qps on different 
regions of the Fermi surface is calculated and the subse- 
quent Raman response is obtained. Section IV compares 
the results of the calculations in Section III to data on 
overdoped Bi2Sr2CaCu208-i-5 (Bi-2212). Section V in- 
vestigates the role of a self consistent treatment of the 
interaction and the role of vertex corrections. Finally 
Section VI summarizes our results and states our conclu- 
sions. 



II. FORMALISM. 

The intensity of inelastically scattered light can be 
written in terms of a differential photon scattering cross 
section as 
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Here — jrnf? is the Thompson radius, wi, cjg are 
the frequency of the incoming and scattered photon, re- 
spectively, and we have set U = kB = 1. 5*^7 is the gen- 
eralized structure function, which is related to the imag- 
inary part to the Raman response function x-yj through 
the fluctuation-dissipation theorem, the second part of 
Eq. (1). n(ijj) is the Bose-Einstein distribution func- 
tion. The Raman response measures "effective density" 
fluctuations, 
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with Tj- the time-ordering operator and the imaginary 
part is obtained by analytic continuation, iuj —^ lo + iO^ . 

The Raman intensity can be represented as a scattering 
off an effective charge density |n5| 



(3) 
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where a is the spin index, 7(k; w/, wg) is the Raman scat- 
tering amplitude, and w/,s is the frequency of the inci- 
dent, scattered photon, respectively. The effective charge 



density is a nonconserving quantity in contrast to the real 
charge . Here the Raman vertices 7 are related to the 
incident and scattering photon polarization vectors e^''^, 
resulting from a coupling of both the charge current and 
the charge density to the vector potential. In general 
the Raman vertex depends non-trivially on both the in- 
cident and scattered photon frequencies. According to 
Abrikosov and Genkin , if the energy of the incident 
and scattered frequencies loi^ujs are negligible compared 
to the energy band gap, 7 can be expressed in terms of 
the curvature of the bands and the incident and scattered 
photon polarization vectors e^^'^ as 



lim 7(k;w/,ws) 
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where terms of the order of 1 — los/loi are dropped. This 
expression is valid if the incoming laser light cannot ex- 
cite direct band-band transitions. However one might 
question the appropriateness of this approximation for 
the cuprates given that typical incoming laser frequen- 
cies are on the order of 2 eV. 

An alternative approach is based on the experimental 
observation that the electronic continua in the metallic 
normal state for a wide range of cuprate materials depend 
only mildly on the incoming laser frequency. Since the 
polarization orientations transform as various elements 
of the point group of the crystal, one can use symmetry 
to classify the scattering amplitude, viz. 
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where $L(k) are either Brillouin zone (BZH, orthogonal 
over the entire Brillouin zone) or Fermi surface (FSH, 
orthogonal on the Fermi surface only) harmonics which 
transform according to point group transformations of 
the crystal . Representing the magnitude but not the 
k-dependence of scattering, the prefactors can be approx- 
imated to be frequency independent and taken as model 
constants to fit absolute intensities. 

Thus we have simplified the many-band problem in 
terms of symmetry components which can be related to 
charge degrees of freedom over portions of the BZ. While 
sacrificing information pertaining to overall intensities, 
we have gained the ability to probe and compare excita- 
tions on different regions of the BZ based solely on sym- 
metry classifications. This can be illustrated by consid- 
ering the various experimentally accessible polarization 
orientations. 

Using an x, y coordinate system locked to the Cu02 
planes, incident and scattered light polarizations aligned 
along x + y,x — y for example transform according to Big 
symmetry, and thus 



$Bj^(k) = cosik^a) - cos{kya) + 



(6) 



where . . . are higher order BZH. Likewise, e-^^'^ aligned 
along x,y transforms as B2g'. 
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^B2g (k) = sm{kxa) sm{kya) + (7) 

The Aig basis function is 

^Aig (k) =00 + 02 [cos{k^a) + cos{kya)] (8) 
+04 cos(fca;a) cos(fcj,a) + a6[cos(2fc2.a) + cos(2fcj,a)] + . . . , 

with the expansion parameters determined via a fitting 
procedure with experiment. The Aig response is not di- 
rectly accessible from experiments and must be obtained 
by subtracting several combinations of the response for 
various polarization orientations. 

By considering the k— dependence of the basis func- 
tions, it is clear that the Big part of the spectra probes 
light scattering events along the k^ or ky axes, B2g probes 
the diagonals, and Aig is a weighted average over the en- 
tire Brillouin zone. In this manner information about the 
momentum dependence of the qp scattering rate can be 
obtained. Since we are interested in probing anisotropy 
effects we will focus on Big and B2g and not consider the 
Aig channel. 

Since the momentum transfer by the photon is neg- 
ligible compared to the momentum of the electrons in 
metals, we are interested in the q = Raman response. 
With the approximations above, the gauge invariant Ra- 
man response is of the form 
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with the renormalized vertex obeying a Bethe-Salpeter 
equation 7 

7(k, iuj, iu! + irt) = 7(k) + ^ X! X! ^^'^ ~ ~ 
xG(k', iw')G(k', iu' + iw)7(k', iu', iu' + iil). (10) 

Here V is the effective qp interaction in the particle-hole 
channel to be discussed below, and N is the number of 
sites. In the absence of vertex corrections to the Raman 
vertex 7(k), the Raman response is given by 

X;\(q = 0,l]) = |5]7'(k) / ^[/H-/(c. + f7)] 



xG"{k,uj)G"{k,uj + n). 



(11) 



The Green's function is given in terms of the self energy 
S by Dyson's equation, G{k,uj) — [w — e(k) — I](k, 
We use a 2-D tight binding band structure £(k) = 
—2t[cos{kxa) + cos(fcya)] -I- At' cos{kxa) cos(fcya) — /x, and 
from now on we set the lattice constant a — I. 



III. PERTURBATION THEORY AND ROLE OF 
SCATTERING. 



A. Other models. 

Before we introduce the NAFL portion of the calcu- 
lation, we start by considering the simple form for the 
Raman response given in Eq. (pj]). We first consider 
non-interacting electrons. In that case, the imaginary 
part of the Green's functions are simply 6 functions, and 
thus the g — > Raman response is given essentially by 
the Lindhard function weighted by the Raman vertex. 
Since for g = there is no phase space available to cre- 
ate electron- hole pairs at low energies, this limit would 
yield no Raman cross section. Thereby, it is crucial to 
have a non-zero self energy in order observe inelastic light 
scattering at all in this limit. 

Zawadowski and Cardona considered free electrons in 
the presence of a potential due to impurities [^. Us- 
ing Eq. (|^ & |l^) and assuming mainly isotropic impu- 
rity scattering, I](k, w) = —i/tl^o, as well as a separa- 
ble interaction diagonal in Fermi surface harmonics (j)L, 
V{k-k') = Vk,k' - EL^Li0L(k)(/.i(k'), the resulting 
cross section can be written in a Lorentzian form: 



1 + i^riy 



(12) 



where a flat band with density of states per spin at the 
Fermi level Np has been assumed, and l/r£ = 1/tl=o — 
1 /tl is the impurity scattering rate in channel L reduced 
by vertex corrections. Finally 7^ is the prefactor of the 
Raman vertex in channel L (see Eq. ^) . While impurities 
do yield a channel dependent response if the Vll are 
different for different L, since the spectra have a peak at 
the scattering rate and then fall off at large oj, it cannot fit 
the large flat continuum observed in Raman experiments 
extending up to the scale of a few eVs . 

A form similar to Eq. ( |l2[) has been obtained for the 
case of qp scattering in a nested Fermi liquid [Q. The 
calculated spectra have the same form as ( p^ ) except that 
a frequency and temperature dependent scattering rate is 
introduced which is of the form 1/t£ = max{(3i,T, aiw}, 
with channel dependent prefactors /3, a. This form for the 
ualitatively good agreement with the ob- 
2^, but does not yield appreciably large 
differences between response functions for different chan- 
nels. In principle features of the qp interaction due to 
nesting, which is pronounced for particular momentum 
transfers near the nesting wavevector, may well account 
for the observed scattering anisotropy. However, calcula- 
tions have not been performed for model Fermi surfaces. 



B. NAFL Theory. 

Therefore, we now consider a form for the qp interac- 
tion due to an effective spin-fermion model interaction 

in 



response gives c 
served spectra 



1 + (q - Q?e - i^l^sf ' 
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This interaction forms the basis of the NAFL theory, and 
has been derived by considering the interaction of elec- 
trons to external spin degrees of freedom [^,0 and de- 
cribes at low temperatures and several dopings the effec- 
tive interaction calculated via FLEX and Quantum 
Monte Carlo treatments of the 2D Hubbard model. 
Recently the optical conductivity was explored within 
this model and was found to be in relatively good agree- 
ment with experiments . Our efforts complement this 
study and investigate the consequences of Eq. ([l^) in 
greater detail. 

The self energy at lowest order is given by 

E(k, ^uJ) = -^J2T.^(^^- P' - *^')G(p, iu;'), 

iuj' p 

(14) 

with G the single particle Green's function. Here LOgf 
and ^ are the phenomenological temperature dependent 
spin fluctuation energy scale and the correlation length, 
respectively, which have been determined via fits to mag- 
netic response data [ [1^ . For optimally doped and under- 
doped systems, the strong AF correlations produce three 
magnetic regions in the normal state. In the phase di- 
agram constructed in Ref. uJsf and ^ obey certain 
relations depending on temperature and doping regimes. 
For each scaling regime, uJsfS,^ = constant (i.e. temper- 
ature independent), with z the dynamical critical expo- 
nent. For the z = 2 regime the spin correlations are 
governed by temperature driven fluctuations at a tem- 
perature above Tcr- As the temperature is reduced be- 
low Tcr the z = 1 or pseudo-scaling regime is reached 
where the spin correlations are strong enough to lead to 
changes from the classical mean field z = 2 regime. Tcr is 
taken to be doping dependent and rises with underdop- 
ing. Therefore the z = 1 region plays a more dominant 
role in the phase diagram than z = 2 for underdoped sys- 
tems while for appreciably overdoped systems, the z = 1 
region vanishes. At lower temperatures still and away 
from overdoped systems a crossover to a pseudogap state 
occurs at T*. However, we will not concentrate on this 
region and defer a discussion of pseudogap effects to the 
last section. 

As in Ref. we use the following approximations: 
(i) we neglect the real part of the self energy, (ii) evalu- 
ate the Green's function at lowest order in the coupling 
constant g^, (iii) neglect all vertex corrections for the 
Raman vertex, and (iv) momentum sums are replaced 
by Ek = /d(^^k/ I Vk I) / de{k). While (iv) does not 
crucially affect the results, approximations (i-iii) - while 
simplifying the calculations - considerably miss impor- 
tant qp renormalizations at larger values of the coupling. 
Therefore we expect that these calculations would be 
most appropriate to describe the spectra taken on over- 
doped cuprate superconductors. 
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FIG. 1. Quasiparticle scattering rate around a section of 
the Fermi surface for different temperatures in the NAFL 
model. Here we have taken ujsf /2t = 0.00876+6.14T/2t for all 
T and ujsf^fit = 0.1058 for T/2t < 0.0273 (z = 1), and then 
taken LOsfi^ /2t = 0.0635 for 0.6 > T/2t > 0.0273 (z = 2). 

1. Quasiparticle scattering rate. 

We first start by consider the effective qp scattering 
rate around the Fermi surface. In Fig. |l| we show the 
qp scattering rate, defined as the imaginary part of the 
self energy S(k, w = 0,T), around the Fermi surface as 
a function of temperature T. Here the effective coupling 
constant is taken as ag^ /2t = 6.25. As a result of our per- 
turbative approach, the magnitude of the scattering rate 
scales directly with ag^ . The band structure parameters 
are chosen to be most applicable to slightly overdoped 
systems: t' /t = 0.45 and a filling < n >= 0.8. The other 
parameters of the theory are chosen to provide a smooth 
link from the z = 2 to the z — I scaling regimes, as 
defined in the figure caption. 

Fig. |l| shows that the scattering rate is extremely 
anisotropic around the Fermi surface, having "hot re- 
gions" (near the zone axes) where the qp scattering rate 
is peaked, and "cold regions" (near the zone diagonals) 
where the rate is smallest. The "hot spots" are those 
regions of the Fermi surface which can be connected by 
Q. These regions are close to the zone axes for the given 
band structure, and the Big geometry most effectively 
probes these hot spots while the i?2g geometry probes 
along the zone diagonals and therefore sees "colder" qps. 
The difference between the scattering at the "hot" and 
"cold" regions is smallest at high temperatures, and is 
due largely to the enhanced density of states near the 
zone axes. As the temperature cools the "hot spots" be- 
come more developed and more strongly peaked as the 
interaction becomes sharply peaked for momentum trans- 
fers near Q. 

Next we consider the frequency dependence of the scat- 
tering rate away from the Fermi level. This is shown in 
Fig. 1^, where the qp scattering rate at T = is plotted 
for both "hot" and "cold" qps, as defined in the figure 
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FIG. 2. Quasiparticle scattering rate for "hot" and "cold" 
spots on the Fermi surface as a function of frequency for 
T = 0. The "hot" spot corresponds to kr = (3.026, 0.4817) 
while the "cold" spot corresponds to kp = (1.107, 1.107). 
Here the z = \ parameters as given in the caption of in Fig. 
|l| are used. 

caption. Both scattering rates vary as uP' at small w, 
characteristic of a Fermi liquid, but cross over to a be- 
havior which is effectively linear in to at larger frequen- 
cies. For the "hot" qps this crossover happens at smaller 
values of u) (0.025, in units of 2i) than for the "cold" qps 
(0.092). Moreover, the scattering is once again larger for 
the "hot" spots than for the "cold" . Similar results have 
been obtained in Ref. [ p^|JT^ , and the reader is directed to 
that reference for further details. From the point of view 
of Raman scattering, we would expect that the differ- 
ent magnitude and frequency/temperature dependence 
of the scattering rate around the Fermi surface will be 
manifest in the channel dependent Raman response. 



2. Channel dependent Raman response. 

Considering the Raman response, we see that the Big 
part of the spectra probes the "hot regions" of large qp 
scattering while i?2g probes the "cold regions". This 
means that the Big response should probe less coher- 
ent qps than the B2g response, which should be reflected 
in the temperature dependence and general lincshapc of 
the Raman continuum. 

Within our perturbative approach we present in Fig. 
^ the results for the Big and B2g Raman response func- 
tions for both z = 1 and z = 2 scaling regimes. The 
parameters we have used for both scaling regimes are 
g = leV, a ^ 3.1 states/eV, t = 250meV, t'/t = 0.45, 
filling < n 0.8, 7Big(k) = bi[cos{k^) ~ cos(fcy)], 
7^2^ (k) — b2 sin{kx) sm{ky) . In addition, for the z — 1 
scaling regime we have used ujsf^ = 50meV and ^"^ = 
0.1 -I- 4.64T/2f, while for the z — 2 scaling regime, 
cj^^^2 ^ eOmeV and ujsf/2t = 0.0237 + 0.55r/2i. These 
parameters are similar to those used in Ref. flj] to de- 
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FIG. 3. Electronic Raman response for the Big 

and Big channels evaluated at different temperatures 
{T/2t = 0.01, 0.02, 0.06 from top to bottom) for the z = 1 
and the z — 2 scaling regimes of NAFL theory. Here we have 
set 6i = 62 = 1, as defined in the text. 



scribe the Hall conductivity data in YBa2Cu307. Here 
the absolute magnitude for the scattering is arbitrary and 
determined by the dimensionless coefficients 61 and 62, 
which are set by fitting to the data. This has no effect 
on the frequency dependent lineshapes. 

The spectra for both scaling regimes share several fea- 
tures. First, the flat continuum at high frequencies which 
is present in the Raman data from all cuprate super- 
conductors is reproduced by the theory. This is a con- 
sequence of a scattering rate S"(ij-') which is effectively 
linearly dependent on lo at frequency scales larger than 
u!sf, as shown in Fig. ^j. The spectra rise linearly with uj 
with a smaller slope for large scattering than for weaker 
scattering. More importantly, the Raman response is dif- 
ferent for the two scattering geometries as a consequence 
of the strong anisotropy of the scattering rate. This is 
due to the relative magnitude and frequency dependence 
of the scattering rate. As shown in Figs. |l| and ||, the 
scattering rate probed in the Big channel is larger than 
B2g for all frequencies and temperatures and becomes ef- 
fectively linear in frequency at a much smaller frequency 
than that probed by the B2g channel. Thus the position 
of the peak of the spectra in the Big channel is at a larger 
frequency in all figures compared to the B2g channel, and 
the Big spectra is overall flatter than the i?2g spectra. 

The differences between the z = 1 and z — 2 results 
are largely quantitative but once again show the influence 
of the anisotropy of the qp scattering. The Big spec- 
tra show the greater difference between the two scaling 
regimes, as the overall spectra appears flatter for z = 1 
than for z = 2. This is due to the more pronounced hot 
spots in the z — 1 regime. The Big channel thus is more 
sensitive than the B2g channel to the growth of the cor- 
relation length as temperature is lowered. Features of a 
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FIG. 4. The Big and -629 Raman continua are plotted 
for T/2t = 0.03 for values of the effective coupling constant 
Qsff = 1, 5, 10, as defined in the text. The z = 1 parameters 
were used for Wsf and ^ as in Fig. ^. 



crossover from z = 2 to z = 1 scaling as the temperature 
is lowered through T^r would be most evident in the Big 
channel pst - 

The role of the effective coupling constant is to de- 
termine the overall shape of the continuum due to the 
increase in the qp scattering rate. An increase of the 
coupling constant leads to (1) an increase in the overall 
magnitude of the cross section, (2) flatter spectra, and 
(3) the movement of the weak peak of the s^pectra out to 
larger energy scales. This is shown in Fig. |j which plots 
the Big and B2g response as a function of frequency shift 
at low temperatures T/2t — 0.03 for different values of 
9eff = ag^/2t. 

As a consequence, the sensitivity of the Raman spectra 
to the details of the parameters can be a very useful tool 
to probe the anisotropic qp scattering rates in much the 
same way as it has been used to probe the anisotropy of 
the energy gap A(k) in the superconducting state p|. 



IV. FIT TO OVERDOPED BI-2212. 

As we remarked, the approximations used (i — iv listed 
in Section III B) are most appropriate for systems with 
weak spin fluctuation scattering and therefore we expect 
that our results would best represent the data from ap- 
preciably over-doped cuprate superconductors. Before 
we consider the effects of a more consistent treatment of 
the interactions in NAFL beyond perturbation theory, it 
is useful to explore how the perturbative theory compares 
to the data on overdoped cuprates in their normal state. 

For a comparison to these data we thus consider the 
results from the z — 2 scaling regime in closer detail. 
Over-doped samples typically have higher residual resis- 
tivities than optimally doped samples. While this may 
be due in part to sample preparation, it is believed that 
the over-doped cuprates increased effective dimension- 
ality allows the qps to interact more strongly with de- 
fects residing out of the Cu02 planes Therefore in 
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FIG. 5. Fit of the 2 = 2 NAFL theory to the Big 
and -B29 spectra taken in Q on appreciably over-doped 
Bi2Sr2CaCu208+a (T^ = 55K) at 60K, 90K, 150K, and 200;t" 
from top to bottom. Parameters used are given in the text. 



addition we will consider an isotropic impurity interac- 



tion H. 



imp 



'k')-Ri J 



Ck' cr, where 



denotes the position of the impurity labeled by i and 
U is the impurity potential. After averaging over the 
position of the impurities, this adds a momentum inde- 
pendent term to the imaginary part of the self energy 
^imp — T^riiNp I C/ 1^, where rii is the impurity concen- 
tration and Np is the density of states per spin at the 
Fermi level. 

The fit of the theory to the temperature dependent 
spectra for each channel obtained in appreciably over- 
doped Bi 2212 (Tc = bhK) by Hackl et al. in Ref. § is 
shown in Fig. for the Big and i?2g channels. The im- 
purity scattering rate for this material can be estimated 
from the extrapolated T = dc resistivity of approxi- 
mately 30/xr2— cm Assuming a Drude resistivity and 
using a plasma frequency of 1.6 eV given in yields 
^imp = 40 cm~^. We note that cugf and ^ have so far not 
been determined via fits to magnetic response data. The 
fits are obtained by choosing: tOsf and ^ at a fixed tem- 
perature, a magnitude of the coupling constant ag^, and 
lastly bi and 62. Then, only the temperature dependent 
part of LUgf was modified to fit the data at other tempera- 
tures. We have used 62/^1 = 0.417, Ws/^^ = 13meV, and 
LUsf = 160iir -I- 0.06T[A']. This parameter choice leaks 
to remarkably good agreement with the data. The flat- 
ness of the continuum for both channels emerges from 
the theory and the relative peak positions come out cor- 
rectly. More importantly, for the first time the tempera- 
ture dependence of the spectra in both channels can be 
explained. NAFL provides an excellent description of the 
qp dynamics in the over-doped region. 



6 



Further information on qp dynamics can be obtained 
from the slope of the spectra at vanishing frequencies. 
We denote this slope as T\ — limQ^o fl / X\{^) , for 
A = Big or B2g. In Ref. Q it was shown that the in- 
verse of this slope has qualitatively different behaviors 
for different doping regimes of various cuprate materi- 
als. While the inverse slope observed in the B2g channel 
always was found to track the temperature dependence 
of the DC resistivity for all regions of doping, the Big 
inverse slope showed a remarkable sensitivity to dop- 
ing. For overdoped systems, the Big slope was similar 
to the i?2g rate, and effectively followed the behavior 



cuprates, Fs^ 



For optimally doped 
constant. Moreover, 



constant -I- 
- T, while Tb,^ 
for underdoped materials the Big inverse slope was found 
to increase with decreasing temperature, indicative of in- 
sulating behavior. 

Within the level of our approximations, the low fre- 
quency Raman response is given by 



(15) 



df^i, 

I Vk 



■7^(k) / 



da; 



[2 cosh(2;/2)]2 j:"{k,xT)' 



The inverse slope samples the qp scattering lifetime at 
regions of the Fermi surface selected by the light polar- 
ization orientations. It is important to note that in this 
level of approximation the momentum dependence of the 
Raman vertices competes with the momentum dependent 
scattering rate on the Fermi surface since the largest con- 
tribution to the integral in Eq. (15) comes from where 



the Raman vertex is the largest and where the qp scat- 
tering rate is the smallest. Therefore, the Raman inverse 
slope washes out direct information about the k— depen- 
dence of the scattering rate and leads to a more similar 
behavior of the Big and i?2g inverse slope, which hides 
the underlying anisotropy of the qp scattering rate. This 
is due to the approximations of restricting consideration 
to the Fermi surface only and to the neglect of vertex 
corrections. In the following section this is shown to be 
crucial if the spin fluctuations are strong. 

However, for weak spin fluctuations, we believe that 
the data on overdoped materials can be best fit with this 
perturbative approach. In Fig. ^ we plot the inverse 
Raman slope obtained from the parameters used in Fig. 
H and compare the results to the data taken in |^ . The 
agreement confirms the fits in Fig. (^ . The data are con- 
sistent with a scattering rate which is not too anisotropic, 
and is governed by elastic scattering from impurities and 
Fermi liquid inelastic scattering from spin fluctuations. 
We remark that the fits are robust to small NAFL pa- 
rameter changes. 
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FIG. 6. Fit of the inverse slope of the Raman response 
for the two symmetry channels. The error bars are from Ref. 
0. 



V. SELF CONSISTENCY AND VERTEX 
CORRECTIONS. 

While good agreement between the theory and the 
data has been shown for over-doped Bi 2212, several fea- 
tures remain to be explained in optimally and under- 
doped systems. As mentioned in the last section, the 
difference of the inverse slope of the Raman response 
between Big and B2g channels grows remarkably upon 
lesser dopings Moreover, features which have been 
associated with a pseudo-gap have been observed in the 
under-doped systems Since our fits to the data 

were restricted by our approximations to weak spin fluc- 
tuations, we need an improved formalism for an insight 
into the effects of strong spin ffuctuations. 

Therefore for the remainder of the letter we solve 
the model lifting the restrictions {i ~ iv) listed before. 
The calculations were performed by flrst solving for the 
Green's function Eq. (|lj) self consistently while main- 
taining a filling (n) = 0.8. We used a hard frequency 
cutoff oj* = 0.4eV for the interaction Eq. (1) to improve 
convergence, but note that our results hardly were af- 
fected by this cutoff. Then the Bethe-Salpeter equation 
( pO| ) was evaluated for the renormalized vertex. The cal- 
culations were carried out on the imaginary frequency 
axis using 256 Matsubara frequencies and 64 x 64 to 
128 X 128 k-points. These values are similar to that used 
for resistivity studies [p6| . The integral equation for the 
vertex 7(k, iuj, icu + ifl) was solved for each external (in- 
ternal) bosonic (fermionic) Matsubara frequency ifl {iuj). 
Once conipleted, the Raman response was evaluated us- 
ing Eq. H and the results were analytically continued 
to the real frequency axis using Fade approximants . 
Since the Raman response is in general quite feature- 
less the Fade procedure was highly accurate (using 20-60 
points changes the results by less than 5 %). 
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FIG. 7. Comparison of the z = 2 (over-doped) results ob- 
tained from perturbation theory (sohd hue), self consistent 
treatment (short dashed line), and the additional inclusion of 
vertex corrections (long dashed lines). 



We first consider how the perturbation results for the 
z — 1 regime are affected. The results are summarized in 
Fig. 1^ which compares the effects of self consistency and 
vertex renornializations to the results from perturbation 
theory for T = O.OSt. ujgf,^, and geff were chosen as in 
Fig. ^ By noting that the position of the peak in the 
Big channel decreases in frequency while the i?2g peak 
position increases under self consistency, it can be seen 
that the effect of a self consistent treatment is to decrease 
the relative anisotropy in the scattering rate around the 
Fermi surface leading to a more effectively "isotropic" re- 
sponse which places the peak of the spectra at nearly the 
same position (~ 0.2t) for both channels. The Big spec- 
trum is slightly narrower in profile than the perturbation 
results, while the B2g spectrum is wider. Although we 
have not tried to fit the self consistent theory to the data 
on over doped Bi 2212 in Fig. |5[ it appears that if dis- 
order is once again added, an equally good fit would be 
obtained using a slightly larger value of the coupling con- 
stant Qef f to offset the reduction due to self consistency. 
We note however that the z — 2 self consistent treatment 
would not yield a particularly good fit to the data on op- 
timally and underdoped systems, which show much more 
anisotropy than that indicated in Fig. ^. 

This can be partially offset by interaction renormal- 
ization of the Raman vertex. When the vertex correc- 
tions are included, the Big spectrum is most affected, 
with an overall reduction intensity occurring over a wide 
region of frequencies, which we interpret as effects re- 
lated to precursor SDW formation . This can be seen 
more clearly by examining the Raman vertices directly. 
In Fig. ^ we plot the Big and i?2g vertices for selected 
directions in the Brillouin Zone for different values of the 
coupling geff and different scaling regimes. The z — 2 
curves correspond to the parameters used in Figs. ^ and 
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FIG. 8. Renormalization of the Raman vertices for the two 
symmetry channels evaluated for the lowest bosonic Mat- 
subara frequency. The solid lines are the bare vertices, 
the dashed are the renormalized vertices for the z = 2 
(over-doped) parameters as in Fig. ^ and the dotted lines 
are the renormalized vertices corresponding to z — 1 (un- 
der-doped) parameters, using a = 16.7 states/eV, corre- 
sponding to YBa2Cu306.5 |l3| . 

0, while the z — 1 curves correspond to the values used 
for YBa2Cu306.5 in |13|, with a = 16.7 states/eV. Since 
g = 1 eV and t = 250 meV has been used for all curves, 
g^fj: = a^^/2t = 6.2(33.4) for the z ^ 2{z ^ 1) curves 



in Fig. ^ We found no cases of symmetry mixing and 
the renormalized vertices retain their bare group trans- 
formation properties. Compared to the bare vertices, the 
Big vertex shows much more of reduction than the i?2g 
vertex. The renormalized vertex 7 for the Big channel is 
reduced in particular near the zone axes or "hot spots" . 
However, the i?2g vertex is hardly affected by the vertex 
renormalizations. This is consistent with the effects of a 
"cold" interaction probed along the zone diagonals. For 
z = 2 the effects of the vertex corrections are small but 
increase dramatically for z = 1. As the interaction is in- 
creased, the Big vertex is reduced dramatically near the 



1 the Big vertex is reduced by a 



"hot spots". For 
factor 2 near the "hot spots" compared to its bare value, 
while the i?2g vertex is reduced by a smaller value 25%, 
and has a peak which is shifted out towards the "hot 
spots" . Since the Raman vertex renormalization enters 
with the opposite sign as the self energy vertex correc- 
tion, these results complement the results of Ref. ]2^ ] 
which showed that the effect of vertex corrections on the 
interaction yields an effectively increased coupling con- 
stant. 

The vertex renormalization suppress the Big cross sec- 
tion compared to the B2g cross section for larger inter- 
action strengths. Since the B2g vertex is affected less 
dramatically, the net effect would appear that the i?2g 
cross section remains relatively unchanged with increas- 
ing interaction, while the Big overall cross section comes 
down in magnitude to be closer to the B2g value. We 
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FIG. 9. The full Raman response for z = 1 (under-doped) 
(evaluated with self consistency and vertex renormalization) 
at temperatures T/2t = 0.03 (sohd), 0.04 (dotted), 0.5 
(dashed), 0.6 (long-dashed), 0.7 (short-dashed), for the Big 
and channels. The parameters used corresponding to val- 
ues for YBazCugOe.s jll. 

remark that this trend is consistent with the data from 
the cuprates which show a large decrease of the Big to 
B2g ratio of the Raman spectra upon under-doping |^] 
which we relate to the growth of "hot spots" intensity 
and subsequently large vertex renormalizations. 

This can be seen more clearly in Fig. which shows 
the full Raman response for larger values of the coupling 
geff corresponding to YBa2Cu306.5 |13|. The spectra 
are much flatter than those calculated for smaller cou- 
pling constants, and appear to differ only quantitatively. 
The effect of the anisotropy of the scattering rate ap- 
pears to have been washed out due to self consistency, as 
the spectra have a "shoulder" at the same frequency and 
show roughly the same temperature dependence. How- 
ever, the overall magnitude of the spectra at larger fre- 
quency shifts is reduced for the Big spectra relative to the 
B2g spectra, leading to a ratio of the continuum nearly 
equal to 1. As indicated in Fig. ||, this is entirely due 
to the vertex renormalization. While the relative mag- 
nitude of the continua and overall lineshape is similar to 
that seen in experiments at a single temperature |^], the 
temperature dependence seen in experiments is qualita- 
tively different. 

This can be explored in more detail by plotting in Fig. 
|l0| the inverse slope of the Raman spectra for the self con- 
sistent and vertex corrected calculation. These results 
are qualitatively different from the perturbative treat- 
ment as discussed earlier. From Eq. (|l^) the inverse 
Raman slope has the same behavior as the qp scattering 
rate at zero frequency. If no disorder is included then the 
inverse Raman slope varies as below a characteristic 
T* and crosses over to T at larger temperatures. T* is 
connected with cosf, and is generally dependent on pa- 
rameter choices. However Fig. |l^ shows that the vertex 
corrected, self consistent treatment does not show sim- 
ilar behavior. For the z—2 scaling regime, the inverse 
slope is proportional to T down to the lowest tempera- 
tures reachable for our level of numerical accuracy, while 
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FIG. 10. The temperature dependence of the Raman in- 
verse slope calculated self consistently with vertex corrections 
evaluated for both symmetries and scaling regimes. 



for the z — 1 regime, the inverse slope appears to ex- 
trapolate to a non-zero value at T = 0. This is also seen 
in FLEX calculations [|lT|. A simple connection between 
these and the perturbative results is that the self con- 
sistency leads to a strong downward renormalization of 
ajsf and a much lower crossover T* to Fermi liquid-like 
behavior. The renormalization is stronger for the z = 1 
case than z = 2, in agreement with our earlier discus- 
sion. Our results thus indicate that the inverse slope 
should fall as T is decreased until ~ Tc- which separates 
the z = 1 and z = 2 scaling behavior, at which point the 
slope should become less T— dependent. Such behavior 
has been reported by Naeini et al. in slightly overdoped 
Lai_2:Sr2:Cu04 and agrees well with the NAFL re- 
sults. 

While this behavior is similar to what has been ob- 
served in the B2g channel, the theory does not predict 
a considerable difference for the Big channel. From our 
calculations, it appears that the effect of self consistency 
obscures the anisotropy inherent in the model, leading 
to a temperature dependent behavior which is similar for 
both channels. While the effect of vertex corrections ap- 
pears to heighten the anisotropy, as shown in Fig. H, 
the self consistency appears to dominate the response, as 
seen in Fig. 
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It would appear then that a treatment 
which neglected self consistency and focused only on ver- 
tex renormalizations would lead to a more apt description 
of the experiments. This remains to be explored. 
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VI. SUMMARY AND CONCLUSIONS 

In summary we have seen how anisotropic qp scatter- 
ing leads to manifestly channel dependent Raman cross 
sections and how through a careful examination of the 
temperature dependence of the spectra, important infor- 
mation can be extracted concerning anisotropic qp dy- 
namics. The coherence or incoherence of the qp's spectral 
function can be probed directly via Raman scattering, 
thereby offering an additional avenue besides ARPES 
to explore momentum resolved qp information. Signa- 
tures of a strong interaction are: (1) the position of a 
weak peak in the spectrum occurring at large frequency 
shifts, (2) the increased overall magnitude of the spec- 
tra, and (3) anomalous temperature dependence of the 
spectra. If the interaction is strongly anisotropic, these 
features are different for different polarization channels. 
The qp dynamics can be probed by examining the tem- 
perature dependence of both the low and high frequency 
portions of the channel dependent Raman response. For 
the cuprates, the consequence of an interaction which 
scatters qps with relative momentum of the antiferro- 
magnetic lattice vector Q are (1) a stronger overall sig- 
nal, (2) a peak at higher frequencies, and (3) a weaker 
temperature dependence of the slope of the response in 
the Big channel compared to the channel. 

Excellent agreement of the NAFL model was found 
with the Raman results on overdoped Bi 2212 (and La 
214 1^^). With the inclusion of impurity scattering, 
the Raman response calculated in a perturbative way 
in NAFL yields a correct description of the spectra at 
a variety of temperatures. Improvements to the theory 
can be made with more accurate determination of the 
parameters of the theory via other experimental probes. 

The effect of a more self consistent treatment leads 
to an overall reduction of anisotropy effects and a more 
isotropic Raman response, in contradiction to the exper- 
iments on optimally- and under-doped cuprates, at least 
for the Big channel. The self consistency yields a much 
weaker Fermi liquid-like behavior than that obtained via 
a perturbative approach. However, the theory does not 
adequately describe the larger relatively incoherent scat- 
tering probed by the Big light orientations. Improve- 
ments can be made by once again exploring different 
parameter choices, but the nearly temperature indepen- 
dence of the Big response seen in the cuprates must be 
obtained by tuning both LUgf and ^ to have a strongly in- 
creasing interaction peaked near (tt, tt) transfers to offset 
the loss of phase space for qp scattering as the tempera- 
ture is lowered. This remains to be explored. 

Interestingly, the role of vertex corrections may help 
to obtain such a behavior. For larger couplings, the ver- 
tex corrections were found to be particularly large near 
the "hot spots" and led to a subsequently reduced vertex 
and overall signal for Big polarizations. Only a small 
change was seen for the B2g channel. This general be- 
havior has been observed in Y 123 and La 214, and has 



been attributed to the opening of a pseudo-gap. More 
than likely it would be important to also employ ver- 
tex corrections to the self energy and to go beyond the 
simple calculations as performed here. Diagrammatically 
we have evaluated the lowest order diagrams for the self 
energy, and we include the interactions in the Raman re- 
sponse through uncrossed vertex corrections (ladder ap- 
proximation). We do not expect that these calculations 
can describe very strong spin fluctuations, but our ap- 
proach in principle should give us an insight as to what 
is the most important physics we might expect as the 
strength of the spin fluctuations is increased by remov- 
ing oxygen. This remains an open challenge. 
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